Energy focusing inside a dynamical cavity 
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We study the exact classical solutions for a real scalar field inside a cavity with a wall whose 
motion is self-consistently determined by the pressure of the field itself. We find that, regardless of 
the system parameters, the long-time solution always becomes nonadiabatic and the field's energy 
concentrates into narrow peaks, which we explain by means of a simple mechanical system. We 
point out implications for the quantized theory. 



The dynamics of confined cavity fields interacting with 
the cavity wall is of great interest for the understand- 
ing of a variety of problems such as hadron bag models 
101, sonoluminescence cavity QED and black hole 
radiations Previous works have mostly approached 
the problem assuming an externally imposed wall mo- 
tion, neglecting the effects of the radiation pressure, or 
used the adiabatic approximation [^j6| . In this paper we 
study, without any approximation, the dynamics of a real 
scalar field inside a cavity, the wall of which moves ac- 
cording to the combined force of a static potential V{R) 
and the field pressure. This system bears important re- 
semblances to more complicated ones, such as the Dirac 
and electromagnetic fields, since they can be partially or 
completely cast in the form of a wave equation. Moreover 
the classical solutions should be a good approximation to 
the quantized fields at least in the case of a large number 
of field quanta. As initial condition for the field we al- 
ways consider a normal mode of the static cavity. This is 
in fact a common situation in the study of many physical 
systems. 

We find that in general the system evolves nonadia- 
batically, and the field energy concentrates into narrow 
peaks. This phenomenon can be understood with the 
help of a simple classical mechanical system. 

In the present work we use natural units and hence the 
action S is dimensionless as are the velocities. This sim- 
ply means that, although we are dealing with a classical 
system, for convenience the action is taken in units of h. 
In one space dimension and with the field only inside the 
cavity, the system is defined by the action 

S = dt' I iMi?2 - V{R) + dx\ [<Pl ^cg]^. 

(1) 

Imposing 5S — Q under any variation of the dynamical 
variables that vanishes a,i t' — Q and t' = t we obtain: 
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Notice the dependence on R of the boundary condi- 
tions. If 4>{R) = is imposed, the total energy, which 
is conserved for a static cavity, is no longer constant for 
Ry^O. Eq. H is satisfied by (l>{x,t) with 



(/)(x,t) ^ G{t - x) + G{t + x) 



(5) 



and the positive sign between the two C?'s ensures that 
the first boundary condition of Eqs. |^ is met. Substitut- 
ing Eq. H in the second of Eqs. |^ we obtain: 



G'{t + R{t)) 



1 - R{t) 
1 + R{t) 



G' {t - R{t)) 



(6) 



For prescribed wall motion, G(z) for any z can be 
found by using Eq. |^ and the null line method |^]. It 
is assumed that the cavity is static for t < to with a 
length R{to). This is equivalent to saying that there is 
a static zone z < zq = to + R(to), in which G(z) is an- 
alytically known. One can find the values of G{z > zq) 
outside the static zone by first solving the algebraic 
equation z = tcqv + R(tcqv) for toqv and then finding 
z^ = toqv ~ Ritcqv)- This proccss, which is equivalent 
to constructing a null line connecting the points z and 
z^, can be repeated many times until a point Zg in the 
static zone is reached. The values of G{z) and G{zs) are 
related through Eq. ^. However in the case under study, 
we do not have, in general, a static zone, and we need to 
verify that knowing the initial conditions of the system 
is enough to implement the above method. 

We will show that in order to find (/((x, t + dt) with 
< X < R(t + dt), it is necessary and sufficient to know 
G{z) and G'{z) for t - R{t) <z<t + R{t) and R{t') for 
t < t' < t + dt. That is just what is required in order to 
have a unique solution of the system of two second order 
equations (||) and (||). 

Since (/)(x,t + dt) = G{t + dt - x) + G{t + dt + x), 
we need to find G{z) and G'{z) for t + dt - R{t + dt) < 
z < t + dt + R{t -f dt). Now we have two cases: either 
z < t + R(t) or z> t + R(t). 
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In the first case it is also true that 

z>t + dt- R(t + dt) >t- R{t) 

as long as i? < I7 *-e. in all physical situations, so that 
we already have the solution. 

In the second case we have to solve the equation 
z = tcqv + R{tcqv), as explained previously. We have 

t + R{t) < icqv + R{tcqv) < t + dt + R{t + dt), 

which, with i? > — 1, implies 

t < ioqv < t + dt. 

Having found t^qv we can derive G'{z) from Eq. ^ be- 
cause, with Zoqv = tcqv " -R(icqv), 

t - R{t) < Zoqv < t + R{t) < 1 , 

so that again we have the necessary information to de- 
termine the evolution of the field. G(z) can then be ob- 
tained by the numerical integration of G'{z). Note how- 
ever that while R — 1 still admits a solution for the field, 
R = —1 doesn't, because the boundary condition requires 
G' {t + dt- R{t + dt)) G' (t + 2dt - R(t)) = 0, which 
in general is inconsistent. Evolving backward in time, 
i.e. with dt < 0, the opposite would be true. 

Using the procedure above we have studied the case 
with V{R) = \K{R - i^o)^ solving, step by step, Eq. | 
numerically by a standard finite difference method. 

As initial condition for the field we choose the funda- 
mental mode of the static cavity with Eqs. ^ as the b.c, 
R{to) = i?o, and i?(to) = 0: 

\ 6 = siliLOtQ COSUJX , TT , , 

Si , w= . (7) 

y (Pt ~ UJ COS WIo COS LOX , Rq 

For convenience we define the dimensionless parameters 
a and f3: a = M/lo, P = Vt/uj — yjK/M /uj, and we set 
the amplitude of the initial field to be 1. In the case of 
a wall initially at rest and with a large mass compared 
to the initial energy of the field, we expect the dynamics 
not to depart considerably from the adiabatic one, that 
is, the wall's motion should be well approximated by the 
solution of Eq. ||, with the field's pressure term replaced 
by its static wall counterpart and the solution of Eq. ^ 
by 

TT 

(j){x,t) — s\n uj{t)t COS, uj{t)x w{t) = — — ■ . (8) 

R[t) 

In order to check the reliability of our numerical im- 
plementation of the algorithm, we first considered a large 
mass of the wall [a = IOOO/tt, (3 = l/(107rV2)). We ver- 
ified that the total energy is very well conserved and the 
motion of the wall is well reproduced by the solution of 
Eq. with the static wall solution for the field pressure. 



We then used a smaller mass keeping K constant, 
i.e. a = IOO/tt and (3 = 1/{-k^/^). As shown in Fig. ||, 
both the wall motion and the field energy density become 
nontrivial. An interesting feature is the concentration 
of the energy density, shown in Fig. |l|c. This is con- 
firmed by the plot of the energy density at two instances 
t = 3A9R{to) and t = 697R{to) in Fig. | compared with 
the static cavity solution. The two peaks at t = 697i?(to) 
move in opposite directions, and their widths decrease 
in time. This phenomenon is even more evident with 
a = IO/tt and /3 = l/(7r\/2) (Fig. ||a), showing a com- 
plex distribution of the peak locations and heights. The 
total energy of the system is the same in all cases. 

Even for the case in Fig. ||b (a = IOOO/tt, f3 — 
l/(107rA/2)), for which we observed the adiabatic evo- 
lution lasting for a long time after to, we can still, letting 
the system evolve long enough, observe the squeezing of 
the field energy density in spite of the slow motion of 
the wall. Keeping K constant we found that the time at 
which the focusing of the energy starts increases roughly 
linearly with M. This suggests that, as one takes into ac- 
count the backreaction of the field on the wall motion, the 
long-time dynamics always becomes nonadiabatic. We 
have verified that this remains true also changing the 
boundary conditions so that the field equals zero at the 
boundaries. 

We believe that the origin of this phenomenon lies in 
the mechanism of energy exchange between the wall and 
the field. To explain it we give the following qualita- 
tive argument. Let's consider the interaction between 
the wave inside the cavity and the wall. At some in- 
stance, the peak of the wave will hit the wall, which can 
be moving either outward or inward. In the former case, 
there will be a transfer of energy from the field to the wall, 
and the speed of the wall will increase slightly. The wave- 
fronts following the peak will lose more and more energy 
to the wall, since the wall moves faster with each succes- 
sive collision. As a result the spatial width of the energy 
distribution decreases. When the wall moves inward, the 
wave gains energy from the wall, and the wavefronts fol- 
lowing the peak gain less because the wall moves slower 
with each successive collision. Again the width of the 
waveform decreases. After some time, this effect leads to 
a drastic concentration of energy into narrow peaks. 

Our argument depends only on kinematics and should 
therefore be applicable not only to waves but many other 
systems, such as a set of particles bouncing back and 
forth in a dynamical cavity. For simplicity we consider 
the dynamics of a set of massless non-interacting parti- 
cles, each having momentum and energy pi, \pi\ (c = 1). 
Inside the cavity they move unperturbed at the speed of 
light. If a particle bounces on the static wall, its momen- 
tum changes sign. The movable wall is subjected to a 
harmonic potential V{R) — {R — Rq)^ . The particle 
momentum p" and the wall velocity v" after an inter- 
action, which is assumed to be instantaneous, are easily 
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derived from energy and momentum conservation: 



V 

p" 



p' + M {v' - 



-4p'/M 
v") , 



1 



(9) 



where v' and p' are the wall velocity and particle mo- 
mentum before the collision. The above equations are 
derived assuming that the sign of p" is always opposite 
to the sign ofp' , which is true as long as the speed of the 
waU is less than 1 and 2M(1 - v') > p' {p' > 0). 

We consider first a set of 1000 particles all with the 
same initial momentum pi = 0.01/i?(to) and a wall ini- 
tially at rest with M = 1000/i?(to) and O = l/i?(to)- 
Already after a few interactions with the wall we could 
observe a regular transfer of energy from the last to the 
first particles to hit the wall. In Fig. ^ we show the 
momenta of the particles after a time t = 3221it!(io) as a 
function of their position. For clarity only positive mo- 
menta are plotted. It is remarkable that the first particle 
to hit the wall has gained more than one tenth of the total 
energy of the system. The above is a very special situ- 
ation which however demonstrates the process of energy 
transfer among particles. 

We then extend this simple mechanical model to the 
case of an infinite number of particles labeled with a con- 
tinuous index k, each having position q{k) and momen- 
tum p{k)dk. In this way we can define an energy density: 



(10) 



£{x,t)= dk\p{k,t)\S{q{k,t) ~ x) 



Not surprisingly £{x,t) satisfies the wave equation in- 
side the cavity. We numerically simulate such a system 
choosing 2000 particles. Initially, we put two particles at 
each of the 1000 uniformly separated sites, and the pairs 
have opposite momenta p{k) = ±10(7r^ cos^ 7rq(fc) -t- 1). 
In Fig. 0b we plot £(0, t)R^{t), which is evidently similar 
to Fig. although the details of the evolution depend on 
how the particles or the field interact with the wall. 

After a long time we observe the formation of many 
smaller peaks in the energy density. Further work is 
needed to understand the problem of the t — > oo evo- 
lution of the system. 

For the scalar field an important situation to study is 
when $1 = tt/Rq, i.e., when the wall motion is in reso- 
nance with the field inside the cavity. We have computed 
the solutions of Eqs. || and ^ for various masses of the 
wall. In Fig. || we plot the wall's position and the field en- 
ergy density at a; = vs. time in the case of a = IOOO/tt, 
R{to) = Rq and ^(^o) = 0.1. In this case we choose 
^0 = ^o/2 so that ip — and the initial functions Eq. |^ 
satisfy the boundary condition Eqs. ^ with R{tQ) ^ 0. 
Besides the beats in the wall motion, two features are 
important. One is the fact that the wall continues to re- 
turn to its initial position after a time T = Rq. This is 
different from the case of non-resonant wall parameters 
where the back reaction of the field changes the frequency 



of the wall motion. Another remarkable effect, as a con- 
sequence, is the appearance of narrow peaks typical of 
a resonantly driven wall motion (|J^,^. This indicates 
the possibility of transferring a large amount of energy 
to the field even with an external, non-resonant, driving 
force [^,|l^ . As long as the frequency of the cavity wall is 
n — tt/Rq, it is enough to push the wall at the instances 
marked by the arrows in Fig. ^, and this frequency de- 
pends on the mass of the wall and can be much smaller 
than ri; increasing the mass decreases the frequency of 
energy exchange between wall and field. This fact might 
help to by-pass the experimental difficulty of achieving a 
resonant driving force, i.e. at frequency 17, on a mirror 
in order to produce high frequency photons [D . 

We have verified that for a small mass, a = IO/tt, the 
wall period remains close to T = 2Rq so that the motion 
is still resonant j^. 

In Ref. it has been shown that the method of null 
lines can also be applied to waves inside an oscillating 
spherical cavity for any value of the angular momentum. 
However, when considering a self-consistent wall motion, 
the spherical symmetry is achieved only in the case of s- 
waves, for which the radial (0) and angular parts can be 
separated. Defining ip = rcj), so that ip satisfies the one- 
dimensional wave equation, we can apply the null lines 
method. The boundary condition for (/>, derived from the 
action similarly to Eqs. ^, is: 



which however for ip translates to: 



(11) 



Ri;,{R{t),t)^^l^^^^-Mm,t) . (12) 

If we want </> to be finite at r = then we must re- 
quire ip — at r = 0, which is satisfied by writing 
TP^G{t-r)- G(t + r). Eq. O becomes: 



G' {t + R{t))-i]G {t + R{t)) 



1 + R 



G' {t ~ R{t)) ^ 7]G {t-R{t)) 



(13) 



with T] = 1/R (^1 + iij . An effective way to solve Eq. |l^ 

numerically for G {t + R{t)) is to define z = t + R{t) and 
to approximate rj and 7 with a constant value between 
z and z — dz for a small enough dz. Integrating Eq. |l^ 
between z and z — dz we obtain: 



G(z) = 



G{z - dz) + 



aVdz 



1 

V 



(14) 



which turns out to be more accurate than standard nu- 
merical integration. 

The force of the s-wave field on the wall is = 
27r_R^(i) [0j — (j)'^] . For i?(io) = we set as initial condi- 
tions for the fields: 
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COS ujtp sin Lor 



4>{r-, to) 



(15) 



where oj ~ 4.4934/i?,(fo) is chosen such that ^(r, ip) sat- 
isfies Eq. |l]with R{t) = 0. As in the ID case, we observe 
the formation of high energy density regions, although in 
3D, this process is much slower. In Fig. ^ we plot the 
energy density at r = R{t) vs. time for a = 5/4.4934 and 
(3 = ^^8/4.4934. These values of the parameters produce 
a completely non-adiabatic evolution. For larger M or 
smaller K we have to evolve the system for a much longer 
time in order to observe the formation of high energy 
peaks. However we have verified that imposing (j) — 
at r = R{t) the peaks appear much earlier and the dy- 
namics is very similar to the one-dimensional situation. 
With resonant wall parameters, = 7r/i?(io), the fea- 
tures observed in ID remain in 3D. With the b.c. Eq. |ll| 
it is also possible to have resonances with equal to the 
difference between the frequencies of the n**^ mode and 
the fundamental mode of the cavity. However such an 
is close to mT/R{to) if n is large, and such resonances are 
not easily distinguishable from the geometric ones |^]. 

In summary we have applied the null lines method to 
study the dynamics of a scalar field inside a cavity whose 
wall is subjected to a harmonic force and the pressure 
due to the scalar field. We have found that the long 
time evolution of the system is always non-adiabatic, re- 
gardless of the parameters of the system. In particular 
there is an interval of time when the field develops nar- 
row packets in energy density that bounce back and forth 
inside the cavity, which can be understood by means of 
a simple mechanical analog consisting of a set of mass- 
less particles bouncing inside a one-dimensional box with 
a movable wall. Such a system confirms our hypothesis 
that the wall motion provides a mechanism of energy 
transfer from low to high energy regions. We have veri- 
fied that the focusing of energy is a robust phenomenon, 
being insensitive to the type of potential for the wall and 
the presence of an external driving force. 

For a quantized field previous works have shown 
that in the case of a prescribed slow wall motion no pho- 
ton production is achieved. Our results strongly suggest 
that the back-reaction of the field may change signifi- 
cantly the evolution of the system. In particular the sec- 
ond derivative of the wall position, which is one of the 
quantities that determine the number of quanta can 
be much larger than in the adiabatic case, as it can be 
seen from the slope of R in Fig.|l|. If the initial number of 
fundamental mode quanta is large, the peaks in energy 
density in the classical solution can imply the production 
of several high energy quanta. 

We have also studied the special situation in which the 
wall frequency is equal to the fundamental frequency of 
the static cavity field. Remarkably the frequency of the 
wall motion does not change due to the field pressure, 
and thus narrow peaks typical of a resonantly driven wall 



motion are produced. A large amount of energy may be 
transferred to the field by providing mechanical energy 
to the wall when the amplitude of the oscillation reaches 
its minimum. This fact might help to by-pass the experi- 
mental difficulty of achieving a resonant driving force on 
a mirror in order to produce high frequency photons M. 



In a further work we would like to address the problem 
of whether periodical solutions are admitted for this kind 
of system and for which values of parameters. 
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FIG. 1. a) Wall position, b) wall velocity, c) energy density 
of the field at a; = for a = IOO/tt, (3 = 1/(77^20). 
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FIG. 2. Spatial distribution of energy density at 
t/R{to) = (dot-dashed), 349 (dashed line), and 697 (solid 
line) for a = IOO/tt, f3 = l/(7rV20). 
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FIG. 4. Classical particles in a dynamical cavity, with 
M = 1000/ R{to), n = l/R{to), and initial momenta 
0.01/R{t(,). a) Particle momenta at t = 3221R{to). b) Gen- 
eralized energy density at x = 0. 




FIG. 5. a) Wall position and b) energy density in a reso- 
nant cavity with a = IOOO/tt and /9 = 1. 
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FIG. 3. Energy density of the field at a; = for: a) 
a = IO/tt, (3 = 1/{ttV2), h) a = IOOO/tt, (3 = l/(107rV2). 
Notice the time intervals. 



OS 



CO 




mto) 



FIG. 6. Energy density in a spherical cavity at r = R{t) 
vs. time for a = 5/4.4934 and (3 = V8/4.4934. 
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